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ABSTRACT 

We study (n, d) and (d, 3 He) reactions on 208 P6 leading to deeply bound pionic 
atoms with one neutron hole being left. We develop a theoretical model to 
calculate cross sections of deeply bound pionic atoms and quasi-elastic pion 
production. The results are compared with data which have been obtained so 
far. We show that the theoretical model describes the data of both reactions 
well. 

1. Introduction 

Since the suggestion of Toki and Yamazaki for the formation of deeply bound 
pionic atoms such as Is and 2p states in heavy nuclei using direct reactions,0l there 
have been a number of experimental attempts to find these states. These states have 
a pionic halo structure due to a strong repulsive pion-nucleus interaction, which saves 
a pion from being absorbed by the nucleus.S'i The (n,p) pion transfer reactions on 
208 Pb at T n = 420MeV were measured by Iwasaki et al. using the Charge Exchange 
facility at TRIUMF.I No evidence for pionic atoms was found, the reason for which 
was interpreted as due to large momentum mismatch at the nuclear surface where the 
pion transfer reaction takes place. A similar conclusion was made for (d, 2 He) reactions 
performed at SATURNE.H It was then concluded that the pion transfer reactions of 
this type were not suited for the formation of deeply bound pionic atoms.! The same 
conclusion was obtained also by Nieves and Oset.i 

Different reactions were studied theoretically by other authors. Nieves and Oset 



studied the reactions induced by leptons such as (7, ir + ) I and (e,e').l They studied 
also (£, A) reactions in order to reduce the momentum mismatch occurring in (n,p) 
reactions by reducing the Q-value due to the S — A mass difference.! Their latest 
suggestion on the use of low energy in-flight pion is also worth-mentioning, in which 
a photon after the direct capture of the in-flight ir~ in deeply bound pionic states can 
be observed.! 

Toki et al., on the other hand, pushed further the use of the neutron and deuteron 
beams, but added another nucleon to participate in the pion transfer process; pion 
transfer followed by one nucleon pick up reactions such as (n, d) and (d, 3 He).0'0 We 
show the schematic diagram for (d, 3 He) reaction leading to pionic bound states in Fig. 
1. These processes have an advantage to be able to choose the incident energy so as to 
make the reaction recoilless. In this kinematics, the pion carries a small momentum 
into the target nuclear system. Although the (n, d) and (d, 3 He) reactions lead to 
complicated nuclear configurations together with the formation of pionic atoms, the 
recoilless condition may make the process experimentally feasible. In fact, the cross 
sections leading to pionic states (n, I) together with a neutron hole states of the same 
orbital angular momentum ('quasi' substitutional configurations, Al = 0) were shown 
to be enhanced.00. 

The first attempt with the use of (n, d) reactions was performed by Trudel et 
al. with the Chargex facility at TRIUMF.Ill The zero degree experimental data was 
reported by Yamazaki at the International Nuclear Physics Conference (INPC)Jl! 
The forward double differential cross section for the 208 Pb(n, d) reaction at T n = 
AOOMeV was measured with the experimental energy resolution of about lMeV. In 
this experiment we observe that the cross section begins to increase gradually from flat 
background at the excitation energy corresponding to about lOMeV below the pion 
production threshold, not at the pion production threshold. This strength below the 
pion production threshold was expected to come from the deeply bound pionic atom 
formation. Matsuoka et al. measured cross sections of {p,pp) reactions to observe the 
deeply bound pionic atoms and obtained the similar spectrum to that of the (n, d) 
reactions.^! Thus, following the theoretical predictions, the experimental search for 
deeply bound pionic atoms by (d, 3 He) reactions was performed in order to improve 



the data quality using the excellent deuteron beam. Yamazaki et al. performed an 
experiment very recently at GSI and reported that they found a clear peak structure 
due to deep, y bo UIld picnic ato m fonnation.0 The a„t between the data and 
the theoretical prediction was found very good and this fact demonstrated that the 
theoretical model describes the reaction very well. 

In this paper we study the one nucleon pick-up reaction leading to deeply bound 
pionic atoms theoretically. The calculated results are compared with data obtained 
recently. In section 2, we describe the formalism of calculating the (n, d) and (d, 3 He) 
cross sections leading to pionic atoms and the pionic continuum states within the 
effective number approach. In section 3, we show general features of the one nu- 
cleon pick-up pion transfer reactions by considering (n, d) reactions. We present the 
calculated results of (n, d) reactions and compare them with the data. The (d, 3 He) 
reaction is studied in section 4. Section 5 is devoted for summary and discussions. 

2. Theoretical Model 

In this section, we describe our theoretical model which was developed to calculate 
the (n, d) and (d, 3 He) reactions leading to deeply bound pionic atoms. Since we are 
interested in the energy spectra only around the pion production threshold, we pos- 
tulate that the spectrum depends on the energy only through three pionic processes; 
(1) bound 7r~ production, (2) quasi-elastic 7r _ production, and (3) quasi-elastic 7r° 
production. There is a neutron hole for ir~ production processes and a proton hole 
for a 7T° production process in the final state. Other processes are assumed to be 
independent of the energy in this narrow region. In the previous publications, we 
have developed the effective number approach to calculate the bound n~ production 
process. 

00 In this section we describe our model including quasi-elastic processes. 
The (n, d) cross sections for the nuclear target are written in the effective number 
approach as 

with 




The first term in Eq. (1) corresponds to the formation of pion bound states and 
the second one to unbound states. The AE in the first term is defined as AE = 
Q + — B w + S n — 3.519MeV for the (n, d) reactions with the pion binding energy 
(B n ), the neutron separation energy (S n ), and the reaction Q- value so as to make a 
peak structure at a resonance energy of each configuration of pion-particle neutron- 
hole state.0 For the (d, 3 He) reactions we need to use Q.787MeV instead of 3.519MeV 
that appeared above.0 The T denotes the width of the bound pionic state. The pion 
binding energies and widths are calculated using the realistic optical potential i. To 
calculate neutron separation energies for deep neutron states we use the potential 
parameters determined by Speth et al. 0. We use experimental values for S n when 
they are available. The E n in the second term is the energy of unbound pions in 
the laboratory frame defined as E n = —Q — S n + 3.519MeV for 7r~ production in 
(n, d) reactions and p n is the corresponding pion momentum. The factor ^E n p n is 
due to the phase volume of the unbound pion. The quasi elastic tt° production is 
included in a similar way. The elementary cross section (da/dQ) nn ^dn is derived 
from the pp — > g?7t + cross sections by using the charge conjugation. In the (d, 3 He) 
reactions, we need to use (da / dQ) d P ^tir+ as the elementary cross section. We show 
the elementary cross sections in Fig. 2.0 

The <pi n and ipj n in N e ff are the pion and the neutron hole wave functions with a 
resultant angular momentum J. The pion wave function is calculated by solving the 
Klein-Gordon equation with the realistic potential S. We use the harmonic oscillator 
wave function for the neutron wave function for simplicity. The spin wave function 
£i ms with averaging over m s is to take care of the possible spin directions of the 
neutrons in the target. The x's denote the initial and the final distorted waves of the 
projectile and the ejectile, which are treated in the eikonal approximation. 

The treatment of deeply bound pionic states is described in detail in previous 
publication. EE We describe here the quasi-elastic pion production part. The scat- 
tering pionic states are obtained by solving the Klein-Gordon equation with the finite 
Coulomb potential and the optical potential V op t, which is written as 



2cuV opt = -47r[6(r) + B(r)] + 4ttV • {L(r)[c(r) + C(r)]}V 



4tt{^-^ v2c W + ^V^^M}, (3) 



where 



6(r) = pi[6op(r) - e w 6i(p B (r) - p P (r))], (4) 



6d = 6d-^(6g + 26?), (5) 



c(r) = pj x [cbp(r) - £ 7r ci(p„(r) - p P (r))], (6) 
B(r) =p 2 B p 2 (r), (7) 



C(r)=p 2 - 1 C p 2 (r) 

47T 

The kinematical factors are defined as 



L(r) = {l + Y X[c(r) + C(r)}}- 1 . (9) 



p^l + w/Af, (10) 

p 2 = l+uj/2M } (11) 
and the reduced energy Co is given by 

u = w/{l+w/AM), (12) 

where u is the pion total energy, e n is the pion charge, A is the nuclear mass number, 
and M is the nucleon mass. The potential is taken from Strieker et al. @, which has 



been obtained by fitting to the pion nucleus elastic scattering cross sections. We show 
the potential parameters in Table 1. Needless to say, but the Coulomb potential is 
strongly attractive for the negative pions, while it is not present for the neutral pion. 
The Klein-Gordon equation is solved by multipole expansion for all the necessary 
partial waves and the relevant pion energies. 

3. Numerical Results on (n, d) reactions 

In this section we will show general features of the one nucleon pick-up pion trans- 
fer reactions by considering (n, d) cases. We would like to first study the threshold 
behavior of the (n, d) reactions leading to 7r~ and tt° states. First of all, we show 
in Fig. 3 the pion wave functions for I = at T v = 5MeV obtained by solving the 
Klein-Gordon equation with and without the optical potential as an example. It is 
important to note that the negative pion feels a strong Coulomb attraction, while 
the neutral pion does not. This difference makes a large difference in the behaviors 
of the pion wave functions in the nuclear region as depicted in Fig. 3. Without the 
optical potential, the neutral pion wave function is the spherical Bessel function (jo). 
The negative pion wave function is pulled in to the small r region relative to the 7r° 
wave function as seen in Fig. 3 (a). By introduction of the optical potential, the 
central components of the wave functions are largely pushed outwards for both the 
negative and neutral pions as seen in Fig. 3 (b). The large difference remains in the 
nuclear interior region. We comment here that these wave functions provide good 
elastic scattering cross sections.il 

These wave functions are then used in eq. (2) for N e ff at various pion energies. 
In order to calculate the energy spectra in the pion unbound region, we have to 
calculate cross sections for many final sates of the final nucleus 207 Pb. As a typical 
case, we show in Fig. 4 the quasi-elastic pion production cross sections for the cases of 
relatively large cross sections for 7r~ and tt° starting from the pion threshold energy. 
Note that the hole states are neutrons for tt~ and protons for 7r°. This figure clearly 
shows the large difference between the threshold behaviors of negative and neutral 
pion formation cross sections. The cross section for negative pion production at the 
threshold is finite, which smoothly connects to the finite cross sections in the bound 



pion region. We make similar calculations for all the configurations of nucleon hole 
and pion particle states until the cross sections converge in the region of interest. The 
single particle states included in the calculations are I/7/2 — 3pi/2 neutron-states (16 
states) and I/7/2 — 3si/ 2 proton-hole states (10 states). 

We show the calculated results in Fig. 5 leading to bound and unbound states. In 
this figure, we assume the energy resolution to be lMeV FWHM which is an actual 
value of the (n,d) experiments. We will see, in the next section, more detail structure 
of the energy spectrum with better energy resolution which is realistic for the (d, 3 He) 
reactions. 

The cross sections in the bound pion region are similar to those published before. 
In the unbound region, we find still a significant amount of strength coming from 
bound pionic atom states coupled with neutron holes at higher excitation energy; 
s — d — g shell with hn/2 intruder, which begins at E x ~ 9MeV from the ground 
state of 207 Pb. On the other hand, the difference between the negative pion spectrum 
and the neutral pion spectrum is remarkable as can be seen by comparing the curves 
with 7r~ and tt°. The negative pion spectrum increases rapidly at the threshold, 
whereas the neutral pion spectrum only gradually increases from the threshold, which 
is 2.7MeV below the negative pion threshold. The small jump seen at 8MeV above 
the threshold in negative pion spectrum is caused by the beginning of the next neutron 
major shell with quasi-elastic pions. Adding all the contributions, we get the total 
cross sections as depicted by the thick solid curve. It is very interesting to note that 
the strength in the bound pion region is almost exclusively due to the bound pionic 
states. The peak appearing at ~ 5MeV mainly corresponds to the configurations 
[1-K ® j^ 1 ] — [2p ® P3/2] and [2p <8> p~J~ 2 ] states. 

Finally we overlay the calculated results on the experimental data, which is done 
in Fig. 6. To do so, we multiply a normalization factor N = 1.7 to the calculated 
results. This factor seems a reasonable value as the uncertainty of the effective number 
approach. We assume the flat background to be 0.28[mb/ sr /MeV] which is not due 
to the pionic contribution. The agreement is remarkable. A small peak structure seen 
in experiment in the bound pion region might be identified to the similar structures 



seen in the calculation. 

4. Numerical Results on (d, 3 He) reactions 

In this section, we compare our theoretical results with the data of the (d, 3 He) 
reactions taken at GSI IllH. We also show the spectrum with better energy resolution 
to predict what should be expected when experiments were performed with better 
energy resolution, and describe the underlying structures of the (d, 3 He) excitation 
functions. We show the energy dependence of the excitation function of the (d, 3 He) 
reaction with the 208 Pb target, too. For (d, 3 He) reactions, we only need to change 
reaction kinematics, elementary cross section, and distortion cross sections used in the 
eikonal approximation for the projectile and the ejectile from (n, d) cases described 
in section 2. 

We show the calculated results on 208 Pb(d, 3 He) reactions at Td = 600MeV in Fig. 
7. In the upper part of Fig. 7, we show the forward cross sections with the resolution 
of 500 keV in order to compare with the experimental spectra shown in the middle 
part of Fig. 7. The agreement is almost perfect. We assume 20[fib/sr/MeV] flat 
background in the theoretical spectra. In order to see the underlying structures and 
also the spectra with a better experimental resolution, we show in the lower part of 
Fig. 7, the cross sections with 200 keV resolution. We can see now many peaks below 
the pion production threshold, which is denoted by the vertical dashed line. The peak 
seen experimentally consists of two states; [l n ® j' 1 ] = [p®Pi/ 2 ] an d \p^P^M- Since 4 
neutrons are allowed in the p^/2 state and 2 neutrons in the py2 state, the ratio of the 
cross sections leading to p^h and p^h is 2. Those two states provide the peak and the 
shoulder structure at Q = —136MeV (around 5MeV binding energy). The pionic s 
state seems to be only weakly excited due to the necessity of some momentum transfer. 
In addition, there is an appreciable amount of contribution from the neutron / 5 / 2 hole 
coupled to the pionic s state to this small peak. Many more peaks are seen closer to 
the pion production threshold. The dominant contributions are again coming from 
the pionic atom states coupled with the p 3 / 2 and pi/ 2 neutron holes. It is interesting to 
comment on the bump structure seen in the calculated spectra around Q = —143MeV 



in the continuum. This bump structure is caused by exciting [ls n £g> s7, 2 ] state, where 
a S1/2 neutron in the 50 < N < 82 major shell is picked up and a pion is placed in the 
Is orbit. In the calculation, the neutron separation energy of the Si/ 2 neutron hole 
state is calculated using the phenomenological potential @ and width of the state is 
assumed to be zero since there is no experimental information. 

We study now the energy dependence of the (d, 3 He) cross sections. As could be 
guessed from the energy dependence of the elementary differential cross sections and 
that of momentum transfer, the (d, 3 He) cross sections at different energies decrease. 
It is, however, interesting to see how the cross sections will change with the inci- 
dent energy. For this purpose, we show the calculated results at Td = 500, 800, and 
lOOOMeV" with the resolution of AE = 200keV and 500keV in Fig.8. As shown in 
the figures, the energy resolution of the cross section is essentially important to see 
the structure of deeply bound pionic states. 

When we lower incident energy to Td = 500MeV, the peak structure at E ~ 
5MeV is enhanced relative to other structure due to smaller momentum transfer. The 
peak structures are almost exclusively produced by the p neutron hole contribution. 
Note that the Is pionic state is not appreciably excited at this energy. Those energy 
dependence are clearly explained by using the energy dependence of the calculated 
effective numbers. As the energy is increased to Td = 800MeV , the cross sections 
are lowered. As indicated by the three kinds of curves in the spectrum with AE = 
200keV, now the contributions from the p neutron hole state become small, while 
the contributions from the i 13 / 2 and /s/ 2 neutron hole states increase. Pionic states 
coupled with / 5 / 2 neutron hole state are excited appreciably. Further increased to 
T d = lOOOMeV^, the peak structures are completely dominated by the zi 3 / 2 neutron 
hole state coupled with various pionic atom states. 

5. Conclusion 

In this paper we have studied the (n, d) and (d, 3 He) reactions on 208 Pb for the 
formation of deeply bound pionic states. We have calculated the formation cross 
sections for various cases and compared them with existing data. Our model is found 
to describe the reaction very well and has good predictive power as proven by the 



fantastic agreement of newly obtained data to our prediction calculated before the 
experiments. We believe the model is a clue for the systematic studies of the structure 
and formation of deeply bound pionic states. 

In order to develop the physics of deeply bound pionic atoms, we need to study 
the structure of pionic atoms more precisely by taking into account the configuration 
mixing effect of pion-particle neutron-hole state, the effect of nuclear deformationil, 
and other effects. They may be important in deriving the pion mass in nuclear medium 
from the experiment^. We also need to study the formation cross sections to other 
target nuclei to look for proper cases for observing the certain states of deeply bound 
pionic statesil. It is also interesting to study the formation of the deeply bound pionic 
states on unstable nucleilzl using the same reaction model. We also mention that this 
direct reaction can be applied to the formation of other mesons such as 77 and u inside 
nucleus to extract their properties at finite density^. 
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Table 1 



Pion-nucleus optical potential parameters determined by Strieker et a/S 



b [fm] -0.046 
61 [fm] -0.134 

B [fm A ] 0.007 + 0.19a 
c [fm 3 ] 0.66 
ci[/m 3 ] 0.428 

Co[fm 6 } 0.287 + 0.93i 

A 1.4 
kplfm- 1 ] 1.4 



Figure Caption 

Fig. 1. Diagram for (d, 3 He) reactions to form pionic bound states with a neutron- hole 
state. 

Fig. 2. (a) The differential cross section (da/dfl) lab at degree in the laboratory 
frame for p(p, d)n + as a function of the incident energy T p , derived from the ex- 
perimental values of (dcr/dQ) c ' m ' for p(p, 7r + )<i.0 The data are taken from references 
The solid line is the phenomenological fit to the data.0 (b) The differential 
cross section (da/dQ) lab at degree for p(d, t)n + as a function of the incident energy 
Td/A[MeV/nucleon], derived from the experimental values for d(p,7i + )t. 00 The 
solid curve is obtained by multiplying 1.7 to the theoretical result of Fearing 10 to 
match the experimental cross sections. 

Fig. 3. The negative and neutral pion wave functions relative to 208 Pb at T n = 5MeV 
for l n = 0. (a) The wave functions without the optical potential, (b) The wave func- 
tions calculated with optical potential, which reproduce the elastic pion scattering 
data. 

Fig. 4. Quasi-elastic pion production cross sections in (n, d) reactions for a negative 
pion with p 3 / 2 neutron hole and those for a neutral pion with a d^/2 proton hole. The 
zero is the threshold energy for each case. 

Fig. 5. 208 Pb reaction cross sections at zero degree leading to deeply bound pionic 
atoms and quasi-elastic pionic states calculated within the effective number approach. 
The horizontal axis, tt Binding Energy, is defined to be OMeV at the 7r~ production 
threshold with the ground state of 207 Pb and it corresponds to Q = — 143.4 [MeV] for 
(n, d) reactions. We take lMeV FWHM as an instrumental resolution. 

Fig. 6. 208 Pb(n, d) reaction cross sections at zero degree at T n = AOOMeV. The solid 
curve is the calculated result with a normalization factor of 1.7. The flat background 
is assumed to be 280\pb/sr/MeV]. We take lMeV FWHM as an instrumental reso- 



lution. The horizontal axis is defined as in Fig. 5. 



Fig. 7. Forward cross sections of the 208 P6(d, 3 He) reaction at T d = 600MeV; (a) 
Theoretical results with the experimental resolution of 500keV, (b) Observed spec- 
trum by Yamazaki et al. 0, and (c) Theoretical results with 200fceV r experimental 
resolution. In (a) and (c) the flat background is assumed to be 2Q[[ib/sr/MeV]. The 
vertical dashed line denotes the ir~ emission threshold energy. 

Fig. 8. Calculated cross sections of the 208 P6(d, 3 He) reaction at (a) T d = 500MeV, 
(b) T d = 800MeV, and (c) T d = lOOOMel/ with 200keV (left) and 500keV (right) 
experimental energy resolution. Solid curves show the full spectra including all con- 
tributions. Dotted curves show the contributions from the neutron 2p states (ps/2 
and P1/2) and dashed curve show those from the neutron i 13 / 2 state. Added is also 
the contribution from the neutron 1/ states (/ 5 / 2 and /V/ 2 ) by dash-dotted curve in 
the middle figure. 
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